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Abstract. The fractal structure of directed percolation clusters, grown at the 
percolation threshold inside parabolic-like systems, is studied in two dimensions via 
Monte— Carlo simulations. With a free surface at y^Cx'' and a dynamical exponent 
z, the surface shape is a relevant perturbation when k<l/z and the fractal dimensions 
of the anisotropic clusters vary continuously with fc. Analytic expressions for these 
variations are obtained using a blob picture approach. 



Since the work of Cardy [1], one knows that a large-scale alteration of the shape 
of a system may influence its local critical behaviour (further references can be found 
in the review [2]). 

Cardy studied systems limited by corners or wedges which are scale invariant 
shapes for an isotropic system and thus, introduce a marginal perturbation leading 
to varying local exponents. More generally, one may consider systems limited by a 
parabolic-like surface at y = ±Cx^ in the (x, j/)-plane in two dimensions, a parabolic 
cylinder or a paraboloid in higher dimensions. Under a uniform rescaling of the lengths 
by a factor 6, C is changed into C — b'^^^C. As a consequence, C grows to infinity for 
fc > 1 whereas it goes to zero when fc < 1 . It remains scale invariant when k = l, which 
corresponds to the marginal corner geometry mentioned above. 

Since 1/C vanishes for the the semi-infinite system, it can be considered as a 
scaling field for the reference fixed point corresponding to the flat surface. This scaling 
variable is an irrelevant perturbation for k>l and a relevant one for k<l. 

In the relevant case, the shape renormalizes to a narrow system, so that the local 
order is weakened. Using a conformal transformation from the half-space to the 
interior of the parabola, one can show that the critical correlation functions decay from 
the tip along the x-axis in a stretched exponential way [3] . The order parameter at the 
tip also displays an essential singularity as a function of the reduced temperature [3, 4, 
5]. Similar results were obtained for the percolation problem [6]. In the case of chain 
polymers (self-avoiding-walks) with one end fixed at the tip, the exponent associated 
with the radius of gyration becomes anisotropic [7] . 
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When the unperturbed system is itself anisotropic, i.e. when the correlation length 
exponents z/y in the a;-direction and in the transverse direction are different with a 
dynamical exponent z given by the ratio /v^ . the appropriate scaling factor becomes 
in the a;-direction. Then the scaling field 1/C transforms as 

— = 6^-^'= - (1) 
C C ^ ' 

and the critical bchavionr is modified when k < 1/z instead of 1. The directed 
walk (polymer) with z = 2 was studied in [8, 9]. In the marginal case, which then 
corresponds to the true parabola, the local susceptibilities given by the numbers of 
walks with N steps starting at the tip and ending (uther anjrwhere or on the surface, 
diverge with a C-dependent power of N as expected. 

In the present work, we study the structure of directed percolation clusters confined 
inside a "parabola". Wc consider the site percolation problem, directed along the 
diagonal of a square lattice, with the sites inside or on the "parabola" belonging to 
the system. 

Using a standard Monte-Carlo method, 2.10^ clusters starting from the tip at 

x~y — Q were generated at the percolation threshold = 0.705489(4) [10] as shown 
in figure 1. The square radii of gyration for s-site clusters along the two directions 

^l = ^.t^-^--^f ys=^.t(y.-y,f (2) 

were put into bins with respect to the sizes 2* < s < 2*+^ [i = 0, 1, 2, . . .) and the 
bin averages (X^), (5^/), were calculated for "finite" clusters, i. e. for clusters with 
x^g^<L = l{){){), the system size along the x-direction. 
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Figure 2. Fractal dimensions rf||(A;), dj_(k) and df(k) versus the system— shape 
exponent k. The hnes give the analytic results obtained via a blob picture approach. 
The fractal dimensions are varying when the perturbation is relevant {k < 1/z = 
0.633(2)). 



Asymptotically, the following power laws are expected 

Sr^X-/'^"^ -Z'"^'^ (3) 

where s is the center of the bin, — {X^y^"^, = {Ygy^-^ and (i||(fc), d^{k), are 
the fractal dimensions in the two directions. In a log-log plot, a linear behaviour is 
obtained for 2® < s< 2^^ 2^''. the deviations at larger sizes being mainly due to the 
finite size of the system. There, the contribution from longer clusters with the same 
values of s is not taken into account. The fractal dimensions, shown in figure 2 versus 
k for different values of C, were obtained from least-square fits of the data in the 
intermediate linear regime. 

In order to analyse our results, let us first consider unconfined directed percolation 
[11] . In two dimensions, the best estimates for the critical exponents, obtained through 
series expansions, are [10]: 

i/j^ = 1.097(2) = 1.734(2) 7 = 2.278(2) (4) 



givmg 



1.581(3) (3 = 0.276(3) (5) 



from scaling. Extending an argument of Stauffer for ordinary percolation [12], one 
obtains the following relations for the fractal dimensions of critical directed percolation 
clusters: 

d = ^±2 = 1.473(2) = = 2.329(3). (6) 
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Figure 3. In the blob picture approach, confined clusters are built through piling up 
successive anisotropic blobs, with a transverse size i„x fixed by the geometry of the 
system. The correlations inside each blob are the same as for unconfined clusters. 



The numerical value oi d^> d may be at first sight surprising for a fractal object. 
Actually, it can be traced to the anisotropy. Let Z|| and Zj_ be the length and width of 

a s-site cluster with s ^ /y'^ii ~ l^"^-^ . It follows that 

^=l^'^^/'^\\ =1^' (7) 

grows faster with s than 1^^. s itself grows faster with Zj^ than an isotropic massive 
object due to the parallel growth. Introducing the characteristic length associated 
with the surface of the cluster 

I ~ (Z||Zx)^^^ ~ s^^'^^ (8) 
a single fractal dimension dj, can be defined with [11] 



and now, = 1.805(2) < 2. 

When the perturbation introduced by the surface is relevant, the structure of 
confined clusters can be studied studied using a blob picture approach [13, 14, 7]. 
Inside the " parabola" , the cluster configuration results from the piling up of anisotropic 
blobs as shown in figure 3. Within each blob, the correlations are the same as for 
unconfined clusters. The nth blob at has a width = Cx^ which is fixed by the 
geometry of the system. Using equation (7), its length is 

Z„ll = C^x^J' ~ C^n- (10) 
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where the last expression corresponds to a continuum approximation in the blob 
index ri. From x,j — x^_i = Z„||, one deduces the following differential equation in 
the continuum limit, 

^ C^ff (j{u-z)/zk^-l+v/zk r^-^^ 

An 

where the last term follows from the n-dependance of in (10). Identifying powers 
of C and n, one obtains u= z/{l — zk) and v = ku. These exponents are positive for 
a relevant perturbation only. Otherwise, when k> 1/z, the cluster grows freely, i. c. 
there is a single blob. Let s„~(^„||)''ii be the number of sites inside the nth blob. For 
the whole cluster, we have: 

N „Ar 
g — y^g ^ C^d.w/O'-zk) I ^zkd\\/(l-zk) ^ (jzd\\/{l-zk) j^[l+zk{d\\-l)]/{l-zk) _ ^-^2) 



In the same way, the cluster length is given by 



Eliminating N between (12) and (13), one finally obtains: 

s~C^('i||-i) d||(fc) = l + 2fc(d|| - 1) Q<k<l/z (14) 

which, together with /j^ = CZ||'^, also gives: 

s~C-i/'= d^{k)='^^ Q<k<l/z. (15) 

The variations of the fractal dimensions are shown in figure 2. For a relevant 
perturbation, (i||(A:) and d_^{k) vary continuously with k until they reach their 
unperturbed values for A; = 1/z. As a result of finite-size effects, the Monte-Carlo 
values deviate from the analytical ones when C increases. For large C values the 
system is more open and the number of blobs inside it remains small for a finite 
system. As a consequence, the Monte-Carlo estimates for the fractal dimensions are 
then closer to the unperturbed ones. 

Results for the tip percolation probability will be discussed elsewhere. 
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